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Abstract — Homology theory has attracted great attention be- 
cause it can provide novel and powerful solutions to address 
coverage problems in wireless sensor networks. They usually use 
an easily computable algebraic object, Rips complex, to detect 
coverage holes. But Rips complex may miss some coverage holes 
in some cases. In this paper, we investigate homology-based 
coverage hole detection for wireless sensor networks on sphere. 
The situations when Rips complex may miss coverage holes are 
first presented. Then we choose the proportion of the area of 
coverage holes missed by Rips complex as a metric to evaluate the 
accuracy of homology-based coverage hole detection approaches. 
Three different cases are considered for the computation of 
accuracy. For each case, closed-form expressions for lower and 
upper bounds of the accuracy are derived. Simulation results are 
well consistent with the analytical lower and upper bounds, with 
maximum differences of 0.5% and 3% respectively. Furthermore, 
it is shown that the radius of sphere has little impact on the 
accuracy if it is much larger than communication and sensing 
radii of each sensor. 

Index Terms — Wireless sensor networks, coverage hole, homol- 
ogy. 

I. Introduction 

WIRELESS sensor networks (WSNs) have attracted con- 
siderable research attention due to their wide potential 
applications such as battlefield surveillance, environmental 
monitoring and intrusion detection. Many of these applications 
require a reliable detection of specified events. Such require- 
ment can be guaranteed only if the target field monitored by 
a WSN contains no coverage holes, that is to say regions of 
the domain not monitored by any sensor. Coverage holes can 
be formed for many reasons, such as random deployment, 
energy depletion or destruction of sensors. Consequently, it 
is essential to detect and localize coverage holes in order to 
ensure the full operability of a WSN. 

There is already extensive work on the coverage hole 
problem in WSNs on 2D plane and 3D space. Some of these 
work used either precise information about sensor locations 
|[T|-|[^ or accurate relative distances between neighbouring 
sensors |[6|-|[8l| to detect coverage holes. The requirement of 
precise location or distance information substantially limits 
their applicability since acquiring such information is either 
expensive or impractical in many settings. Thus connectivity- 
based approaches are of great interest for us. In this category, 
homology-based schemes have received special attention be- 
cause of its powerfulness for coverage problems in WSNs. 
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Homology was first introduced by Christ and his collabo- 
rators in |[9l-fTT| to address the coverage problems in WSNs. 
They introduced a combinatorial object, Cech complex, which 
fully characterizes coverage properties of a WSN (existence 
and locations of holes). Unfortunately, this object is very 
difficult to construct even if the precise information about 
the relative locations of sensors is available. Thus, they in- 
troduced another more easily computable complex, Vietoris- 
Rips complex (we will abbreviate the term to Rips complex 
in this paper). This complex is constructed with the sole 
knowledge of the connectivity graph of the network and gives 
an approximate coverage by simple algebraic calculations. 
More precisely. Rips complex may miss some special coverage 
holes. It is thus of paramount importance to analyse the 
coverage holes missed by Rips complex in order to evaluate 
the accuracy of homology-based coverage hole detection. 

The main contributions of our paper are as follows. First, 
the relationship between Cech complex and Rips complex in 
terms of coverage hole on sphere is analysed. Furthermore, 
it is found that a hole in a Cech complex missed by a Rips 
complex must be surrounded by a spherical triangle. Based on 
that, a formal definition of spherical triangular hole is given. 

Second, we choose the proportion of the area of spherical 
triangular holes as a metric to evaluate the accuracy of 
homology-based coverage hole detection. Such proportion is 
analysed under a homogeneous setting and it is related to the 
communication and sensing radii of each sensor. Three cases 
are considered for the computation of such proportion. For 
each case, closed-form expressions for lower and upper bounds 
of the proportion are derived. 

Third, extensive simulations are performed to evaluate im- 
pacts of communication and sensing radii, radius of sphere 
on proportion of the area of spherical triangular holes. It is 
shown that simulation results are well consistent with the ana- 
lytical lower bound, with a maximum difference of 0.5%, and 
consistent with the analytical upper bound, with a maximum 
difference of 3%. Furthermore, simulation results show that 
the radius of sphere has little impact on the proportion when 
it is much larger than communication and sensing radii. It 
further indicates that our analytical results can be applied to 
more general 3D surfaces. 

The rest of the paper is organised as follows. Section II 
presents the related work. In Section III, the network model 
and the formal definition of spherical triangular hole are given. 
Closed-form lower and upper bounds for proportion of the area 
of spherical triangular holes under three different cases are 
derived in Section IV. Section V compares simulation results 
and analytical bounds. Finally, Section VI concludes the paper. 
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II. Related work 

We present the related work in terms of two aspects: 
coverage ratio analysis and homology-based coverage hole 
detection. 

First, extensive research has been done to analyse coverage 
ratio in a WSN on 2D plane or 3D surface. In fT2|, p3|, 
the fraction of the area covered by sensors was analysed. 
In p4| , the authors studied how the probability of coverage 
changes with the sensing radius or the number of sensors. In 
pSj , a point in a plane is defined to be tri-covered if it lies 
inside a triangle formed by three nodes, and the probability of 
tri-coverage was analysed. In |16|, the authors proposed the 
surface coverage model and derived the expected coverage 
ratio under stochastic deployment on 3D surface. In pTj , 
the expected coverage ratio under stochastic deployment on 
3D rolling terrains was derived. Nevertheless, none of these 
research considered spherical triangular holes. We provided 
some initial results about the proportion of the area of tri- 
angular holes on 2D plane in | [T8| . In this paper, we aim to 
extend results in |T8| on sphere. 

As for homology-based coverage hole detection, De Silva et 
al first proposed a centralized algorithm that detects coverage 
hole via homology in | [T0| . They constructed the Rips complex 
corresponding to the communication graph of the network 
and determined the coverage by verifying whether the first 
homology group of the Rips complex is trivial. Then the above 
ideas were first implemented in a distributed way in |19|. 
It is shown that combinatorial Laplacians are the right tools 
for distributed computation of homology groups and can be 
used for decentralized coverage verification. In f20l, a gossip- 
like decentralized algorithm for computation of homology 
groups was proposed. In | |2T| , a decentralized scheme based 
on Laplacian flows was proposed to compute a generator 
of the first homology group. But all these homology-based 
algorithms do not consider the cases that Rips complex may 
miss some special coverage holes. In this paper, we describe 
such cases and analyse the proportion of the area of coverage 
holes missed by Rips complex. 

III. Models and definitions 

Consider a collection of stationary sensors (also called 
nodes) deployed randomly on a sphere with radius R ac- 
cording to a homogeneous Poisson point process with intensity 
A. For any two points pi and p2 on S^, the distance between 
them d{pi^p2) is defined to be the great circle distance, which 
is the shortest distance between any two points on the surface 
of a sphere measured along a path on the surface of the sphere. 
As usual, isotropic radio propagation is assumed. All sensors 
have the same sensing radius Rs and communication radius 
Rc on It 

means for any sensor located at v on S , any 
point p on with d{v^p) < Rs is within the sensing range 
of the sensor; and for any two sensors located at v^, Vj on S^, 
they can communicate with each other if d{vi^Vj) < Rc. In 
addition, we assume Rg <^ R, Rq <^ R. 

Before defining the two combinatorial objects, known as 
Cech complex and Rips complex, it is necessary to give 
a brief introduction to some tools used in the paper. For 



further readings, see p2}-p4|. Given a set of points V, 
a /^-simplex is an unordered set [vq, ^i, ^fc] ^ V where 
Vi 7^ Vj for all i 7^ j, k is the dimension of this simplex. 
The faces of this ^-simplex consist of all (^-l)-simplex of the 
form [^0, •••7 •••7 ^fe] for < i < /c. For example, 

on a sphere S^, a 0-simplex [vq] is a vertex , a 1-simplex 
[vo,^i] is the shorter arc of the great circle passing through 
^0 and Vl, a 2-simplex [^07^17^2] is a spherical triangle 
VQV1V2 with its interior included, see Figure [T] An abstract 
simplicial complex is a collection of simplices which is closed 
with respect to inclusion of faces. A /c-dimensional abstract 
simplical complex /C is an abstract simplicial complex where 
the largest dimension of any simplex in /C is /c. 



Vo 




0-simplex 1-simplex 2-simplex 

Fig. 1. 0-, 1- and 2-simplex 

Let V denote the set of sensor locations in a WSN on 
with radius R and 5 = {s^,vgV} denote the collection of 
sensing ranges of these sensors: for a location v, Sy = {x ^ 
: d{x,v) < Rs}- Then Cech complex and Rips complex 
can be defined as follows. 

Definition 1 (Cech complex). Given a finite collection of 
sensing ranges {s^, v G V}, the Cech complex of the col- 
lection, C(V), is the abstract simplicial complex whose k- 
simplices correspond to non-empty intersections of k + 1 
distinct elements of {sy^ v G V}. 

Definition 2 (Rips complex). Given a metric space (S^,(i), 
a finite set of points V on and a fixed radius e, the Rips 
complex ofV, lZe(V), is the abstract simplicial complex whose 
k- simplices correspond to unordered (k +l)-tuples of points in 
V which are pairwise within distance e of each other 

According to the definitions, the Cech complex and Rips 
complex of the WSN, respectively denoted by Ci^^(V) 
and 1Zr^{V), can be constructed as follows: a ^-simplex 
[vo.vi,--- ,Vk] belongs to Cr,{V) whenever nfLgS^^ ^ 
and a ^-simplex [vq^vi^ - ■ ■ , Vk] belongs to 1Zr^{V) whenever 
d{vi^Vm) < i^c for all < / < m < k. In addition, since we 
consider only coverage holes on the sphere it is sufficient 
to construct 2-dimensional Cech complex and 2-dimensional 
Rips complex of the WSN, denoted as C^rI{V) and 7^^^^^(V) 
respectively. 

Figure |2] shows a WSN, its Cech complex and two Rips 
complexes for two different values of Rc. Depending on the 
relation of Rc and Rs , the Rips complex and the Cech complex 
may be close or rather different. In this example, for Rc = 
2Rs, the Rips complex sees the hole surrounded by 2, 3, 5, 6 as 
in the Cech complex whereas it is missed in the Rips complex 
for Rc = 2.bRs. At the same time, the true coverage hole 
surrounded by 1,2,6 is missed in both Rips complexes. 

In fact, as proved in |[25|, any coverage hole can be found 
in Cech complex. Furthermore, there are following relations 
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(c) 



(d) 



Fig. 2. (a) a WSN, (b) Cech complex, (c) Rips Complex under Rc 
(d) Rips Complex under Rc = 2.5Rs 



2Rs 



between C%]{V) and 7^'^^^^(V). 

Lemma 1. Let V denote the set of node locations in a WSN 
on with radius R, all nodes have the same sensing radius 
Rs and communication radius Rc, Rs <^ R^ Rc <^ R, then 



(2). 



<^(v)cc£^(v)c7^i^(v), 

whenever Rc < R slt ccos {[3 cos^ {Rs / R) — l]/2) 



(1) 



Proof: See the Appendix. ■ 
According to ([T]), some relationships between Cech complex 
and Rips complex in terms of coverage hole can be derived 
as illustrated in the following corollaries. 

Corollary 1. When Rc < R diTccos{[3 cos'^ {Rs / R) - l]/2), if 



there is no hole in 7^^^(V), there must be no hole in C^r^{V). 

(2) 

Corollary 2. When Rc > 2Rs, if there is a hole in 1V^\V), 

- (2) 

there must be a hole in C^^ (V). 

Corollary 3. When R diTccos{[3 cos'^ {Rs / R) - l]/2) < Rc < 

f2) 

2Rs, there is no guarantee relation between '^^^(V) and 
C^R^(y) in terms of holes. 

- (2) 

From the discussions above, a hole in a C^ (V) not seen 

(2) ^ 

in a (V) must be bounded by a spherical triangle. Based 
on this observation, a formal definition of spherical triangular 
hole is given as follows. 

Definition 3 (Spherical triangular hole). For a pair of com- 
plexes C^^l (V) and 7^^^^ (V), a spherical triangular hole is an 
uncovered region bounded by a spherical triangle (2-simplex) 
which appears in n^^\v) but not in C^r]{V). 



IV. Bounds on proportion of spherical triangular 

HOLES 

In this section, the conditions under which any point on 
with radius R is inside a spherical triangular hole are first 



given. From the discussions in Section [Ill| it is found that the 
proportion of the area of spherical triangular holes is related to 
the relation of Rc and Rs . Three different cases are considered 
for the proportion computation. For each case, closed-form 
expressions for lower and upper bounds of the proportion are 
derived. 

A. Preliminary 

Lemma 2. For any point on S^, it is inside a spherical 
triangular hole if and only if the following two conditions 
are satisfied: 

1) the great circle distance between the point and its closest 
node is larger than Rs. 

2) the point is inside a spherical triangle: the convex hull 
of three nodes with pairwise great circle distance less 
than or equal to Rc. 

Lemma 3. If there exists a point O which is inside a spherical 
triangular hole, then Rs < R arccos ^^[1 + 2 cos{Rc/ R)]/3. 

Proof: According to Definition [3] if there is a point O 
inside a spherical triangular hole, then there exists a 2-simplex 



a e nf^iy) while a ^ e^l{V), so TZf^{V) <t CrI{V). 
According to ([l]), we have Rc > R diTccos{[3 cos^ {Rs / R) — 
l]/2) ^ Rs < arccos y^[l ^ 2cos{Rc/R)]/3. ■ 

Lemma 4. Let O be a point inside a spherical triangular hole 
and I denote the great circle distance between O and its closest 
neighbour, then Rs < I < R arccos + 2 cos(i?c/^)]/3. 

The proof is similar as that of Lemma [T] 

Since we assume nodes are distributed on according to a 
homogeneous Poisson point process with intensity A, any point 
has the same probability to be inside a spherical triangular 
hole. This probability in a homogeneous setting is also equal 
to the proportion of the area of spherical triangular holes. 

We use spherical coordinates {R^O^ip) to denote points 
on with radius R, where is polar angle and (p is 
azimuth angle. We consider the probability of the point N 
with spherical coordinates (i?, 0, 0) being inside a spherical 
triangular hole. Since the communication radius of each sensor 
is at most R^ only the nodes within Rc from the point TV can 
contribute to the spherical triangle which bounds a spherical 
triangular hole containing N. Therefore, we only need to 
consider the Poisson point process constrained on the spherical 
cap C{N^Rc) which is also a homogeneous Poisson process 
with intensity A, where C{N^Rc) denotes the spherical cap 
centered at point N and the maximum great circle distance 
between N and points on the spherical cap is Rc. We denote 
this process as ^. In addition, T(x, z) denotes the property 
that the point N is inside the spherical triangular hole bounded 
by the spherical triangle with points x^y^z as vertices. When 
no,ni,n2 are points of the process ^, T(no,ni,n2) is also 
used to denote the event that the spherical triangle formed 
by the nodes no,ni,n2 bounds a spherical triangular hole 
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containing the point N. In addition, we use T\no^ni^n2) 
to denote the event that the nodes no,ni,n2 can not form 
a spherical triangle which bounds a spherical triangular hole 
containing the point N. 

Let To = ro{^) be the node in the process $ which is 
closest to the point N. There are two cases for the point N to 
be inside a spherical triangular hole. The first case is that the 
node To can contribute to a spherical triangle which bounds a 
spherical triangular hole containing the point N. The second 
case is that the node tq can not contribute to any spherical 
triangle which bounds a spherical triangular hole containing 
the point N but other three nodes can form a spherical triangle 
which bounds a spherical triangular hole containing the point 
N. So the probability that the point TV is inside a spherical 
triangular hole can be defined as 

p{X) = F{N is inside a spherical triangular hole} 

= P{ U T(no,ni,n2)} 

{no,ni ,n2}C<^> 

= P{ y r(To,ni,n2)} (2) 

{ni,n2}C<^)\{ro(<3E))} 

+ P{ U T(no, ni, ns) | (tq, ns, n^)} 

{-0>--->-4} 
C<i.\{To('i')} 

B. Case < Rc < R diTccos{[3 cos'^ {Rs / R) - l]/2) 

Theorem 1. When < Rc < R diTccos{[3 cos'^ {Rs / R) - 
l]/2), p(A) = 0. 

Proof: According to when < Rc < 

means 

that there exists not a spherical triangle which appears in 
n^^\v) but not in C^r\v), so there is no spherical triangular 
holes which means p{X) =0. 

■ 

C Case R8iTccos{[3cos^{Rs/R) - l]/2) < Rc < 2Rs 

Theorem 2. When R diTccos{[3 cos^ {Rs / R) - l]/2) < Rc < 
2Rs, pi{X) < p{X) < Pn(A), where 



Pi{X) =27rX^R^ / sinOodOo / dcpi I smUi 

JRs/R J27V-iPrn ^ Oq 

^ g-A|C(Ar,i?^o)|g-A|5+(^o,^i)|(l _ g-A|5-(^o,^i,^i)|)^5/^ 



(3) 



and 



Pu{X) =27tX'^R^ / sinOodOo / d(pi / sin(9i 

X g-A|C(Ar,i?0o)|g-A|5+(^o,^i)|(l _ g-A|5-(^o,^o,^i)|)^5|^ 

+ P{ U T(no,ni,n2) | (tq, na, 714)} 



C<i>\{xo(<i>)} 



and 



arccos y [1 + 2 cos(i?c/^)]/3 
arccos[(cos(i?c/^) — cos^ 6>o)/(sin^ Oq 
min{l9i^i,6>i^2} 

COs(i?c/^) 



7ini = arccos 



sin^ 6q sin^ 



arctan(cos (/^i tan 6>o) 



0" 

cos{Rc/R)/ sin^ 6>o sin^((/9i - (prn) 



hu2 = arccos 

+ arctan(cos((/:?i — (prn) tan6>o) 
Proof: We first prove the lower bound. It can be obtained 



from ^ that 

p{X) > P{ 



y T(ro,ni,n2)} 

{ni,n2}C$\{ro($)} 



So for the lower bound, we only consider the first case 
that the closest node tq must contribute to a spherical triangle 
which bounds a spherical triangular hole containing the point 
N. 

Using spherical coordinates, we assume the closest node tq 
lies on (i?, ao, 0) and use l^*! to denote the area of the set S, 
then we can get the distribution of ao as 



Fao{0o)=P{ao<0o) 



-X\C{N,R9o)\ 



(5) 



since the event ao > Oq means that the spherical cap 
C{N,R6o) does not contain any nodes from the process, 
which is given by the Poisson probability e-A|C'(Ar,i?6'o)| 
Furthermore, \C{N, R6o)\ can be given as 

\C{N,ROo)\ = r [ R^sinOdifidO = 27rR^{l-cosOo) 
Jo Jo 



From ([5]) and ([6]), we can get the density of tq 

Fc,{dOo) = 2^XR^ smOoe-^\^^^^^^'^\dOo 



(6) 



(7) 



The integration range for Oq can be easily obtained. 
According to Lemma [4} we have Rs < ROq < 

arccos v^[l + 2 cos(i?c/^)]/ 3, so Rs/R < Oq < Oqu = 
arccos ^"[1 + 2 cos(i?c/^)]/3. 

Therefore the probability of the first case can be given as 



P{ U T(ro,ni,n2)} 

{ni,n2}C<3E>\{ro(<3E>)} 



P{ U T((i?,^o,0),ni,n2)}F,oW) 



(8) 



(4) 



where is the restriction of ^ in C{N, Rc)\C{N, ROq). 

Once the node tq is determined, the other two 
nodes must lie in the different half spaces: one in 
= R+ X (0,7r/2) X (7r,27r) and the other in 
H- = R+ X (0,7r/2) X (0,7r). Assume ni lies in 
and n2 lies in H~ . Since the great circle distance to tq 
is at most R^ ni and n2 must also lie in the spherical 
cap C{ro^Rc). Furthermore, the great circle distance to the 
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point N is at most Rc and larger than Ra^, they should 
also lie in the area C {N ^ Rc)\C {N ^ Rao) . Therefore, ni 
must lie in f]C{ro, Rc) i^CiN, Rc)\C{N, Rao) and 
712 must lie in H' f]C{ro, Rc) O^i^^ ^c)\C{N, Rao). In 
addition, considering the great circle distance between 
rii and n2 should be at most Rc and the point N 
should be inside the spherical triangle formed by tq, 
ni and n2, rii must lie in the shadow area = 

n C{ro, Rc) n C{N, Rc)\C{N, Rao) f] C{M2,Rc), 
shown in Figure |3] Mi and M2 are two intersection points 
between bases of spherical caps C{N^Rao) and C{ro^Rc), 
such intersection points must exist in this case since 
Rc < 2Rs < 2Rao. 




Fig. 3. Illustration of area in the case arccos( [3 cos ^(i?s/i?) — 
l]/2) <Rc< 2Rs 

Ordering the nodes in by increasing azimuth angle 
SO that Ti = {R^Oi^ipi) has the smallest azimuth an- 
gle (fi. And assume the nodes tq, n and another node 
r2 e H-f]C{To,Rc)f]C{N,Rc)\C{N,Rao) can form a 
spherical triangle which bounds a spherical triangular hole 
containing the point N, then r2 must lie to the right of the great 
circle passing through ri and N, denoted by H~^{(pi) which 
contains all points with azimuth angle (p G {(pi — ir^ipi). In 
addition, the great circle distance to ri is no larger than R^ so 
the node T2 must lie in the region 5", as illustrated in Figure 

SI 

S- (ro, ri) = S- {ao, O^^i) = H' f] C(ro, Rc) 

f] C{N, Rc)\C{N, ROo) n ^+(^1) n ^(^1 ' 
Here we need to obtain the density of node ri . Considering 
the way n was defined, there should be no nodes with a 
azimuth angle less than (pi in A+, that is to say no nodes 
are in the region 

^+(ro,ri) = = 

Since the intensity measure of the Poisson point process in 
spherical coordinates is XR'^ sin OdOdif, the density of ri 
can be given as 




Fig. 4. Illustration of areas 5'+ and S in the case 

RaTCCos([3cos'^(Rs/R) - l]/2) < Rc < 2Rs 



Fr,{dOi,d(pi) = Ai?^sin(9ie-^l^^("°'^^)ldi9id(/?i (9) 

Then we derive the integration domain D{ao) with 
respect to parameters (Oi^ifi). Assume the point M2 
has the spherical coordinate {R^0o^(pm)^^m ^ (O^tt). 
Since the great circle distance between tq and M2 is 
Rc, then according to the spherical law of consines, we 
have cos{Rc/R) = cos^ ^0 + sin^ Oo cos (pm ^ (fm = 
arccos[(cos(i?c/^) — cos^ 6>o)/(sin^ 6>o)]. It can be seen 
that points Mi and Q have the spherical coordinates 
{R^Oo^^n — (prn) and {R,6o^2(prn) respectively, where Q 
is one intersection point between bases of spherical caps 
C{N^Rao) and C{M2,Rc). Thus the integration range for 
(fi is [27T — iprm2ipm]- In addition, assume any point with 
great circle distance Rc to tq has the spherical coordinate 
{R,Ot,(pt), still using the spherical law of consines, we 
have cos{Rc/ R) = cos ^0 cos Of + sin 60 sin 6t cos (pt =^ 
Ot{(pt) = arccos [cos (i?c/ R)/ V'^ ~ sin^ ^0 sin^ cpt] + 
arctan(cos(/:?^ tan6>o). Similarly, assume any point 
with great circle distance Rc to M2 has the 
spherical coordinate {R^O[^ip[), we can obtain 

9[{(p^) = arccos [cos (i?c/^)/ ^1 ~ sin^ Oo sin^((/pj — (fm)] + 
arctan(cos((/?J — (/:?^) tan 6>o). Then the integration range 
for 01 is [6o,0iu], where 6>i^ = min{6>i^i, 6>i^2}, = 
0t{(pi),0iu2 = 0[{(pi). Furthermore, |5'+(ao,^i)| can be 
expressed as 

\S^{ao,ipi)\= r I R^ sin OdOdcp 

Assume only ro,ri and nodes in 5'~(ro,ri) can contribute 
to the spherical triangle which bounds a spherical triangular 
hole containing the point N, we can get a lower bound of the 
probability that the point N is inside a spherical triangular 
hole. It is a lower bound because it is possible that ri can 
not contribute to a spherical triangle which bounds a spherical 
triangular hole containing point N, but some other nodes with 
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higher azimuth angles in the area can contribute to such 
a spherical triangle. For example, in Figure [5] if there is no 
node in S~ but there are some nodes in S'~ , then ri can not 
contribute to any spherical triangle which bounds a spherical 
triangular hole containing point N, but r[ can form such a 
spherical triangle with tq and another node in . Based on 
the assumption, we have 



Next we will prove the upper bound. As discussed in Section 




Fig. 5. Illustration of areas S + and S in the case 

Rd,Tccos{[3cos'^{Rs/R) - l]/2) < Rc < 2Rs 



P{ U T((i?,^o,0),ni,n2)} 
>P{ U T((i?,^o,0),ri,n2)} 

n2C<3E>',^p5-(ro,ri) 
= //P{ U T((i?,^o,0),(i?, ^1,(^1), 712)} 



(10) 



[f P{*',„(5-(0o,^i,^i)) > OjFrAdOud^i) 

JJd{9o) 

J Jd(0„) 



where \S {9o,0i,ifi)\ can be expressed as 

rf 

cos{Rc/R) — cos^i cos^o 



5'-(i9o,6>i,(/?i)| = / / sin 02d02dif 2 (11) 



and 



^21 = ^1 — arccos ■ 
O2U = min{6>i^ii, ^2n2} 



sm Ui sm 6/0 



i72w2 = arccos 



s{Rc/R)/ ^1 - sin^ 6>o sin^((/92 - ^1) 



+ arctan(cos((/?2 — ^1) tan 6>i ) 

Therefore, from (|7]), ([S]), ([9]) and ([TO]), the lower bound 
shown in ^ can be derived. 



IV-A[ there are two cases for the point N being inside a 
spherical triangular hole. As for the second case that the 
closest node tq can not but some other nodes can contribute to 
a spherical triangle which bounds a spherical triangular hole 
containing the point N, it is not easy to obtain a closed- 
form expression for such probability, but we can get it by 
simulations. Simulation results show that this probability is 
less than 0.16% whenever Rc < 3Rs with any intensity A. So 
we still focus on the probability of the first case. 

For the lower bound, we only considered the case that ri 
contributes to a spherical triangle which bounds a spherical 
triangular hole containing point N. For the upper bound, we 
need to further consider the case that ri can not but some 
other nodes in can contribute to such a spherical triangle, 
shown in Figure |5] Assume the node r{ = {R^O[,(p[) with 
the second smallest azimuth angle in ^4+ can contribute to 
such a spherical triangle, it means that there is no node in 
S~{ao,6i^(pi) but there is at least one node in the region 

Then the density of the pair (ri,r{) is given as 



^^r[{dOi,d(pi,dO[,d(p[) 

= X^R^ sinOi sine'^e-^\^^^'^'^'^'^^\deidifide[dif'^ 



(12) 



The probability that ri can not but r{ can form a spherical 
triangle which bounds a spherical triangular hole containing 
point N with tq and another node in S'~ {aQ^9i^(^i^9[^(^'i) 
can be given as 

P{ U T((i?,^o,0),r(,n4) | ((i?, ^0, 0), n, na)} 

{r.3,r.4}C<i.' 

^0 

p|s-(-o:-;) 



P{^'^o(^"(^o,^i,(^i))=0} 

xp{^'^o(^"(^o,^i,^i,^;,^i))>o} 

,t[ {dOi ,d(pi, dO[ , d(p[ ) 
= JJJ ^e"^l'^"^^°'^i'^i^l X (1 - e"^l'^'"^^°'^i''^^'^i'^i^l) 

Fr,^ri{dOud^udO[,d^[) 

(13) 

As we can see from Figure |5] as long as r[ has a 
higher polar angle than ri has, the sum of \S~ {ao^Oi^(pi)\ 
and \S^~{ao^6i,(pi^6[,(p[)\ will be always smaller than 

\S~{ao,ao,(fi)\. 



Therefore we can get from ( 13 ) 



P{ U T((i?,^o,0),r{,n4) | ((i?, ^0, 0), n, na)} 

P|s-(xo,xp 



< 



-A|5-(0o,^i,^i)| _ ^-A|5-(0o,^o,^i)| 



F,^^,^{dei,d^i,de[,d^[) 



(14) 
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Let 5^+(^o,(^i,(/^;) = 5+(^o,(/^;)\>S+(^o,(/^i), then 



1 _e-A|A+\5+(^o,^i)| ^ I 



(15) 



It is the complement of the probabiHty that no node is in 
the area A^\S^{Oo,(pi). 

From ([To]), ([12]), ([14]) and ( [T5] ), we can obtain 

P{ U T((i?,^o,0),ri,n2)} 

+ P{ U T{{R, ^0, 0), r{ , n4) | ((i^, ^o, 0), n, ns)} 

{ri3,ri4}C<^>' 
^0 

P|S-(xo,r;) 

(16) 



where \S {6o^6o^(pi)\ has the similar expression as (11). 

Similarly, we can further consider the case that neither 
of Ti and r[ can contribute to a spherical triangle which 
bounds a spherical triangular hole containing point N, but 
Other nodes with even higher azimuth angle can contribute to 
such a spherical triangle. In this way, we can get the same 
result as ([16]). 

Therefore, it can be derived that 



P{ U T((i?,^o,0),ni,n2)} 

{ni,n2}C<^>'^^ 

< ff (l-e-^l^"(^°'^°'^^)l)F,,(d^i,d(^i) 



(17) 



From ([2]), ([T]), ([5]), ([12]) and ( [TT] ), the upper bound shown in 
([4]) can be derived. 



D. Case Rc > 2Rs 

Theorem 3. When Rc > 2Rs, pi{X) < p{X) < Pn(A), where 



27rA^i?'^| J^"" sinOodOo J dcpi sin(9i 



X e 



-X\C{N,Reo)\^-X\S^{eo,^i)\^l _ g-A|5-(^o,^i,^i)| 



+ / sin^o<^^o / difi / single" 

^ ^-X\S+ {00,^1)1^1 _ g-A|5-(^o,^i,^i)|)^5|^| 



■X\CiN,Reo)\ 



and 



27tX'^R^^ j^"" sinOodOo J dipi sini9i 



^ g-A|C(Ar,i?^o)|g-A|5+(^o,^i)|(l 

+ / sinOodOo / / 



-A|5-(0o,^o,^i)| 



X e 



'Rc/2R 

-x\s+{eo,^i)\ 



(1 



-A|5-(^o,' 



smt/ie 



\\C{N,R9o)\ 



P{ U T(no,ni,n2) |T'(To,n3,n4)} 



{"O.--- ."4} 
C<1.\{to(<1.)} 



(19) 



and 

^lu = ™Il{^lnl, ^ln2} 

^in2 = arccos cos{Rc/ R)/\/l — sin^ 6>o sin^ {(pi — tt) 

+ arct an (cos ((/:)! — tt) tan^o) 
In this case, we can use the same method as in Section 



(18) 



IV-C to get the lower and upper bounds, shown in (18) 
and ( [T9| ) respectively. But we need to consider two situ- 
ations Rs/R < Oo < Rc/{2R) and Rc/{2R) < Oq < 
Oou = arccos y^[1 + 2 cos(i^c/^)]/3. In the first situation, 
^0 ^ Rc/{2R) means that the spherical cap C{N^R9q) 
is completely included in the spherical cap C{ro^Rc). The 
illustrations for the areas A+, 5+, S~ ^ S'~^ and S'~ are shown 
in Figure |6];a) and[6];b) respectively. In addition, the integration 
range for (pi is [7r,27r]. The second situation is the same as 
that in Section IIV-CI 

V. Simulations and performance evaluation 

In this section, simulation settings are first given. Then 
simulation results are compared with analytical lower and 
upper bounds under different settings of Rg^Rc^ R. 

A. Simulation settings 

A sphere centered at the origin with radius R is considered 
in the simulations. The probability of the point with spherical 
coordinate (i?, 0,0) being inside a spherical triangular hole 
is computed. Sensors are randomly distributed on the sphere 
according to a homogeneous Poisson point process with in- 
tensity A. The sensing radius Rg of each node is set to be 
10 meters and communication radius Rc is chosen from 20 
to 30 meters with interval of 2 meters. Let 7 = Rc/Rg, then 
7 ranges from 2 to 3 with interval of 0.2. In addition, A is 
selected from 0.001 to 0.020 with interval of 0.001. For each 
pair of (A, 7), 10^ simulations are run to check whether the 
point with spherical coordinate (i?, 0, 0) belongs to a spherical 
triangular hole. 

B. Impact of Rg and Rc 

As illustrated in Section |lll| Rg ^ R and Rc <C R, here 
we choose R = lORg to analyse the impact of Rg and Rc on 
the probability of any point being inside a spherical triangular 
hole. Under this configuration, the probability p{X) obtained 
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(a) 




(b) 

Fig. 6. Illustrations of areas in case Rc > 2Rs. (a) the areas A+, 5'+ and 
S- (b) the areas and S'' 



by simulations is presented with the lower and upper bounds 
in Figure [TJa) and [TJb) respectively. 

It can be seen that for any value of 7, p(A) has a maximum 
at a threshold value Ac of the intensity. As a matter of fact, 
for A < Ac, the number of nodes is small. Consequently the 
probability of any point being inside a spherical triangular hole 
is relatively small too. With the increase of A, the connectivity 
between nodes becomes stronger. As a result, the probability 
of any point being inside a spherical triangular hole increases. 
However, when the intensity reaches the threshold value, the 
probability is up to its maximum. p{X) decreases for A > Ac. 
The simulations also show that Ac decreases with the increase 
of 7. 

On the other hand, it can be seen from Figure |7ja) and|7jb) 
that for a fixed intensity A, p{X) increases with the increases 
of 7. That is because when Rg is fixed, the larger Rc is, the 
higher is the probability of each spherical triangle containing 
a coverage hole. 




- simulation y - 2.0 

- lower bound y - 2.0 

- simulation y - 2.2 

- lower bound y - 2.2 

- simulation y - 2.4 

- lower bound y - 2.4 

- simulation y - 2.6 

- lower bound y - 2.6 

- simulation y - 2.8 

- lower bound y - 2.8 

- simulation y - 3.0 

- lower bound y - 3.0 



0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 
intensity k 



(a) 





• simulation y - 2.0 






• upper bound y - 2.0 






w simulation y - 2.2 






w upper bound y - 2.2 






□ simulation y - 2.4 






□ upper bound y - 2.4 






simulation y - 2.6 






upper bound y - 2.6 






simulation y - 2.8 






upper bound y - 2.8 






A simulation y - 3.0 






A upper bound y - 3.0 







0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 
intensity A, 

(b) 

Fig. 7. Proportion of the area of spherical triangular holes under R = 10 Rs 
(a) simulation results and lower bounds ; (b) simulation results and upper 
bounds 



Furthermore, the maximum probability increases quickly 
with 7 ranging from 2.0 to 3.0. These results can also provide 
some insights for planning of WSNs, which will be discussed 
in Section IV-DI 

Finally, it can be found in Figure [TJa) that the probability 
obtained by simulation is very well consistent with the lower 
bound. The maximum difference between them is about 0.5%. 
Figure [TJb) shows that probability obtained by simulation is 
also consistent with the upper bound. The maximum difference 
between them is about 3%. 



C. Impact of R 

Although we assume Rs <^ R and Rc <C R, to better 
understand the impact of R on the probability of any point 
being inside a spherical triangular hole, we choose R to be 
bRsj lORs and lOOi^g. In addition, we also want to know the 
difference of the probability under spherical and 2D planar 
cases. Therefore, simulation results, lower and upper bounds of 
the probability under spheres with radii 5i?s, lOi^s, lOOi^s and 
2D plane are shown in Figure [8];a),[8jb) andjSjc) respectively. 

It can be seen from Figure [8] that simulations results, lower 
and upper bounds under spheres with radii bRgj lORgj lOOi^g 
and 2D plane are very close with each other. More precisely, 
the maximum difference of simulations results under spheres 
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4 5 


6 




X 10"^ - 


— A — sphere R 




— H — sphere R 




— — sphere R 


= 100/?^ 


O 2D 





0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 

intensity A, 



(a) 




0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 
intensity A, 



(b) 




0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 
intensity X 

(C) 

Fig. 8. Comparison of the proportion of the area of spherical triangular 
holes (a) comparison of simulation results; (b) comparison of lower bounds; 
(c) comparison of upper bounds 



with radii bRs and 10i?s is about 0.045%, which is about 
0.06% under spheres with radii bRs and lOOi^s and is about 
0.03% under spheres with radii lOi^g and lOOi^g. In addition, 
the maximum differences of simulation results between 2D 
planar case and spherical cases with radii bRg-, lOi^s, lOOi^s 
are 0.05%, 0.03% and 0.02% respectively. It means the larger 
the radius of sphere is, the more closer are the simulation 
results under sphere and 2D plane, it is because the larger the 
radius of sphere is, the more likely of the local of each node 
on the sphere to be planar. 



With respect to lower and upper bounds, it is found that 
under any two spheres with radii 5i?s, lOi^g, lOOi^g, the 
maximum difference of lower and upper bounds are 0.06% 
and 0.12% respectively. Furthermore, under spheres with radii 
5i?5, 10i?5, 100i?5 and 2D plane, the maximum difference of 
lower bounds is also 0.06%, and that of upper bounds is also 
0.12%. More importantly, under sphere with radius lOOi^g 
and 2D plane, the maximum difference of lower bounds is 
5 X 10-^ and that of upper bounds is 2.5 x 10 ^. It means the 
probabilities under cases of sphere with radius lOOi^g and 2D 
plane are nearly the same, which is quite logical since when 
the radius of sphere is much more larger than the sensing 
radius of any node, the local of any node can be considered 
to be planar. 

It can be further found that under above cases, the maximum 
differences of simulation results, lower and upper bounds 
are all so small that they can be neglected. Consequently, it 
also means that the radius of sphere has little impact on the 
probability of any point on the sphere to be inside a spherical 
triangular hole. Our results can thus be extended to more 
general 3D surfaces. 

D. Discussions on applications 

In this paper, we only consider spherical triangular holes, for 
non-spherical triangular holes, we assume they can be detected 
and covered by additional nodes. Under this assumption, our 
analytical results can be used for planning of WSNs. For 
example, a WSN is used to monitor a mountain and the 
ratio 7 = 2, according to the analytical upper bounds, we 
can see that the maximum proportion of the area of spherical 
triangular holes under 7 = 2 is about 0.06 %, which can be 
neglected. It means that as long as the surface of mountain can 
be spherically triangulated by nodes, we can say the mountain 
is covered. But if 7 = 3 and at least 95% of the surface of the 
mountain should be covered, then it means that the proportion 
of the area of spherical triangular holes can be at most 5%. 
From the analytical upper bounds of 7 = 3, it can be seen 
that when the intensity A = 0.009, the upper bound is about 
5%, SO in order to cover at least 95% of the mountain, the 
intensity of nodes should be larger than 0.009. Similarly, our 
results can also be used in the scenarios when using satellite 
to cover the whole earth. 

VI. Conclusions 

This paper studied the accuracy of homology-based cov- 
erage hole detection for wireless sensor networks on sphere. 
First, the situations when Rips complex may miss coverage 
holes were identified. Then we chose the proportion of the 
area of coverage holes missed by Rips complex as a metric to 
evaluate the accuracy. Three different cases were considered to 
compute the accuracy. For each case, closed-form expressions 
for lower and upper bounds were derived. Simulation results 
are well consistent with the derived lower and upper bounds, 
with maximum differences of 0.5% and 3% respectively. In 
addition, simulation results also show that the radius of sphere 
has little impact on the accuracy as long as it is much larger 
than communication and sensing radii of each sensor. This 
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means that our results can be applied to more general 3D 
surfaces although the results are derived on sphere. This 
problem will be investigated in our future work. 

Appendix 
Proof of Lemma[T] 

Proof: The second inclusion is trivial because for any 
/c-simplex [vq, ^i, • • • , v/c] G C^^^(V), it means the sensing 
ranges of these nodes have a common intersection, so the 
pairwise distance d{vi^Vj) < 2Rs for all < z < j < /c, 
which means ['^q, '^i, • ' ' j^k] ^ (^)- 

^ (2) 

As for the first inclusion, it is clear that 7^}^^(V) and 

- (2) 

^Rs (^) contain the same 0-simplices. It is also easy to see 
that all 1-simplices in T^^^V) must also be in C^^^(V) since 
for any 1 -simplex [vi^vj] with distance d{vi^Vj) < Rc < 
RdiTccos{[3cos'^{Rs/R) - l]/2) < R diTccos{2 cos'^ {Rs / R) - 
1) = 2Rs, it means that the sensing ranges of the two nodes 
have a common intersection. So we only need to prove that all 
2-simplices in 7^^^(V) must be in C^^^(V). It is equivalent to 
say that for any three nodes with pairwise great circle distance 
no larger than Rc, their sensing ranges must have a common 
intersection. 

(2) 

Assume a 2-simplex [^0,^1,^2] ^ ^r^C^)^ then the three 
nodes ^0,^1 and V2 must determine a plane a. We consider 
the spherical cap on cut off by the plane a. Since Rc < R, 
the spherical cap must be on a hemisphere. It is easy to see 
that the intersection of the plane a and sphere is a circle 
c. Let Oi be the center of circle c, O be the center of S^, P 
be the intersection of line OOi and S^. 

Using spherical coordinates, we assume thepoint P has a 
spherical coordinate (i?, 0, 0). P may be insidqjor outside the 
spherical triangle V0V1V2, which is shown in Figure |9ja) and 
[9jb) respectively. 



assume /3 > 27r/3. According to the spherical law of con- 
sines, we have cos(/3) = '"^^°^if2|^^7ff ^"^"^^ < -1/2 ^ 
cos{doi/R) < [3cos'^{dp/R) - l]/2. In addition, doi < 
Rc < RdiTccos{[3cos'^{Rs/R) - l]/2) cos{doi/R) > 
[3cos^{Rs/R) - l]/2, and < doi/R,dp/R < 7r/2, so we 
have [3cos'^{Rs/R) - l]/2 < [3 cos'^ {dp / R) - l]/2 ^ dp < 
Rs, which means the point P is a common intersection of 

- (2) 

sensing ranges of '^o^'^i and V2, so ['^q, '^i, '^2] ^ C^^(V). 

If P is outside the spherical triangle V0V1V2, as shown 
in Figure [9ja), it indicates that the spherical triangle V0V1V2 
must be contained in half of the spherical cap. Assume vq^ vi 
and V2 have spherical coordinates {R,6,(po)^ {R^^^^i) and 
(i?, 6>, (/?2), where G (0, 7r/2), (/?o < ^1 < ^2, then we have 
(fi - if if 2 - ^1,^2 - ^0 ^ (0,7r). Using doi,di2,do2 to 
denote the pairwise great circle distances between vo^Vi,V2, 
then according to the spherical law of consines, we have 



cos{doi/R) = cos^ - 
cos{di2/R) = cos^ - 
cos{do2/R) = cos^ - 



- sin^ cos((/9i — (fo) (20) 

- sin^ ^ cos((/92 — ^1) (21) 

- sin^ 6 cos((/?2 — fo) (22) 



In addition, we use a to denote the angle between two arcs 
and vqV2, M to denote the middle point of the arc vqV2 
and d^Mi <^im to denote great circle distances between vq, vi 
and M. It can be seen d^M = <^02/2. Similarly, we have 



cos{di2/R) - cos{doi/R) cos{do2/R) 

cos a = : — — rn:^ — : — — rr^^ (^-3) 



sm{doi/R) sm{do2/R) 



diM doi doM 

cos —I— = cos — — cos — — 
R R 2R 



. doi doM 
sm — — cos — — cos a (24) 
R 2R 



From (|23l) and ([24|), we can obtain 



/V>. 




cos 



diM _ cos{doi/R) cos{di2/R) 



R 



2cos{do2/{2R)) 



Consequently 



cos 



~r7 



-cos 



~Rr 



cos 



dpi 



cos - 



cos 



(25) 



2cos(do2/(2i?)) 



(26) 



(a) 



(b) 



Fig. 9. Illustrations of P and spherical triangle V0V1V2: (a) P is inside the 
spherical triangle voviV2', (b) P is outside the spherical triangle voviV2. 

It can be seen that P has the same great circle dis- 
tance to ^0,^1 and V2, denoted by dp. If P is inside the 
spherical triangle V0V1V2, as shown in Figure [9ja), then we 
can prove dp < Rs. Since P lying inside the spherical 
triangle V0V1V2 means /3 + 7 + J = 27r, there must be 
one angle no smaller than 27r/3. Without loss of generality, 

^It also includes the case that P is on one arc of the spherical triangle 

V0VIV2. 



From ([20]), ([21]), ([22]) and (|26]), we get 



d 



cos 



IM 



d{ 



R 



- cos 



OM 



sin^ cos ^^^^ sin sin ^^^^ 



R 



cos 



do2 
2R 



(27) 

Since < (fi — (f0i^2 — ^11^2 — ^0 < tt and < 
diM / R^doM / R^do2/ R < 7r/2, it can be obtained from ( [27| ) 
diM < doM < ^c/2 < Rs, which means the point M is a 
common intersection of the sensing ranges of ^o^^i and V2, 



so [^0,^1,^2] ^ (V). It means all 2-simplices in Tl^^^ {V) 

- (2) 

must be in C^^ (V). Consequently the first inclusion is proved. 



(2). 
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